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Heavy ion collisions at RHIC are well described by the (nearly ideal) hydrodynamics. In the 
present paper we study propagation of perturbations induced by moving charges (jets) on top of the 
expanding fireball, using hydrodynamics and (dual) magnetohydrodynamics. Two experimentally 
observed structures, called a "cone" and a "hard ridge", have been discovered in dihadron correlation 
function with large-pt trigger, while "soft ridge" is a similar structure seen without hard trigger. 
All three can be viewed as traces left by a moving charge in matter, on top of overall expansion. 
A puzzle is why those perturbations are apparently rather well preserved at the time of the fireball 
freezeout. We study two possible solutions to it: (i) a "wave-splitting" acoustic option and (ii) a 
"metastable electric flux tubes". In the first case we show that rapidly variable speed of sound 
under certain conditions leads to secondary sound waves, which are at freezeout time closer to the 
original location and have larger intensities than the first wave. In the latter case we rely on (dual) 
magnetohydrodynamics, which also predicts two cones or cylinders of the waves. We also briefly 
discuss metastable electric flux tubes in the near-T c phase and their relation to clustering data. 

PACS numbers: 



I. INTRODUCTION 
A. Outline of the problem 

Experimental data obtained in heavy ion collisions at 
RHIC were shown to be in very good agreement with 
hydrodynamical description of the excited system cre- 
ated, supplemented by hadronic cascades for description 
of hadronic phase and freezeout [TJ |2j |3] . Dissipative ef- 
fects from viscosity provide only small corrections at the 
few percent level [H 16] . In the last few years RHIC ex- 
periments have focused more on two and three-particle 
correlations, which revealed rather rich phenomenology 
of correlations. These correlation phenomena indicate 
that certain fluctuations occuring on top of this overall 
expanding fireball have quite puzzling dynamics, which 
is the subject of this paper. 

Starting the introduction, let me note that issues to 
be discussed are somewhat similar in nature to what 
happened in cosmology in the last decades. While ex- 
istence of background radiation had dramatically con- 
firmed the existence of Big Bang in 1960's, it is only 
more recent observations of the temperature fluctuations 
seen in top of the overall Hubble expansion which made 
cosmology more quantitative science. In fact the angular 
size and magnitude of such fluctuations seen on the sky, 
corresponding to the sound modes at the plasma neu- 
tralization time, now provide excellent data on timing of 
the cosmological expansion. One may hope that studies 
of small perturbations propagating on top of expanding 
"Little Bang" - the heavy ion collisions - will turn out 
to be equally valuable in understanding many details of 
its complicated physics. 

We will turn to observations in the next subsection: 
but before we do so, let me formulate the main physics 
dilemma of this work: either 



(A) the perturbations are indeed small- amplitude hy- 
drodynamical perturbations, yet propagating differently 
from what is naively expected on the basis of a geometric 
optics, 

or (B) they are not hydrodynamical but include certain 
fields/structures, affecting their expansion. 

In this work we will examine whether both of those 
solutions are viable. The option (A) - to be referred 
to as the "acoustic solution" - will reveal creation of 
the secondary waves, induced by time-dependent speed 
of sound. (In fact this effect was already noticed by 
Casalderrey and myself in Ref. [8 in connection with 
conical flow.) As we will show, such secondary waves 
are brighter and smaller in size, as sketched in Fig{2|c). 
However, as we will find, it is not clear whether solution 
A will be viable quantitatively, as it require rather sharp 
drop in a speed of sound. 

The second option B also leads to double cones, now as 
two components of Alfven waves in a (dually) magnetized 
medium. Furthermore, some of them have small or even 
zero expansion velocity, and indication to existence of 
stabilized electric flux tubes in near-Tc temperature in- 
terval. Metastable microscopic flux tubes in the near-T c 
region had also been considered in a different context 
before, by Liao and myself [35l [38] in connection with 
lattice data on lattice potentials and charmonium sur- 
vival. Yet again, although such tubes have good reasons 
to exist, the final conclusion on whether they are robust 
enough to explain the observed "cone" and two "ridges" 
would require a lot of further experimental and theoret- 
ical work. 

Magnetohydrodynamics is a well known part of plasma 
physics, developed by Alfven, Fermi, Chandrasekhar and 
many others since 1940's, see standard textbook such as 
[7]. It is an approximation which keeps only magnetic 
field in Maxwell eqns, while the electric field is assumed 
to be totally screened. Ideal MHD approximation is the 
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limit of infinite conductivity of plasma a — > oo, similar 
to zero viscosity approximation for ideal hydrodynamics. 
In MHD the coupling between the field and and matter is 
obtained by inclusion of the (magnetic) field contribution 
into the stress tensor of the medium. 

Early stages of heavy ion collisions are believed to be 
described by the so called "glasma" , a set of random 
color fields created by color charges of partons of the two 
colliding nuclei at the moment of the collision. As argued 
by McLerran and Venugopalan [11], the color charges 
of partons at particular location in the transverse plane 
cannot be not correlated (they come from different nucle- 
ons) , thus they add up randomly. For large nuclei those 
charges and fields can become large enough to be treated 
classically. However as two discs with charges move away 
from each other, those classical field are getting smaller 
and (in a still poorly understood process) rather quickly 
create the quark-gluon plasma, in which the occupation 
numbers are thermal, O(l), with the classical fields dis- 
appearing. Indeed, perturbative theory of asymptotically 
hot QGP predicts perturbative electric screening mass 
Me ~ gT from the one-loop perturbative polarization 
tensor calculated long time ago [12]. Such perturbative 
approach however creates no screening of the static mag- 
netic fields, Mm = 0, like in the QED plasma. 

Transition from field to plasma has been subject of 
multiple works, too many to be mentioned. One notable 
idea is existence of instabilities [14] such as the so called 
Weibel instability, producing filamentation of the longi- 
tudinal flow and transverse magnetic field B±. Asakawa 
et al [15] argued that such chaotic magnetic field would 
remain in plasma, leading to "abnormal viscosity". Un- 
fortunately those ideas, derived perturbatively and basi- 
cally taken from electromagnetic plasma contexts, cannot 
possibly be valid for QGP at RHIC. 

The crucial difference between the QED and QCD 
plasmas lies in the existence of magnetically charged 
quasiparticles - monopoles and dyons - leading to 
nonzero magnetic screening mass Mm ~ g 2 T first sug- 
gested by Polyakov 30 years ago [13] and by now well 
confirmed by subsequent lattice studies. Thus hot QGP, 
unlike the electromagnetic plasmas, screen both the elec- 
tric and the magnetic fields at some microscopic scales, 
although at a bit different ones. 

Furthermore, recent theory developments known as 
"magnetic scenario" [31] [32] for the near-T c region show 
that the situation with electric and magnetic screening 
masses get in fact inverted in this region. Multiple lat- 
tice studies, e.g. ref [16] (see also Fig. 16 my review [29] , 
where these issues are discussed in more detail) have 
shown that at T < 1.4T C the relation between electric 
and magnetic screening masses get inverted, Mm > Me- 
As T — > T c the electric screening mass strongly decreases, 
partly because of heavy quark and gluon quasiparticles 
and partly because of their suppression by small Polyakov 
loop expectation value < L >, going to zero below T c . As 
temperature decreasing toward T c the magnetic screen- 
ing mass is instead increasing monotonously, to a rather 
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M M (T T c ) * 3T C , M E (T -+ T c ) -+ (1.1) 

Such behavior of screening as well as other theoreti- 
cal considerations had lead to the "magnetic scenario" 
[3Tj [32] . which basically views the near-T c region, from 
about 0.8T C up to 1.4T C as magnetic plasma, dominated 
by (gluomagnetic) monopoles. The near-T c region, for- 
mally known as the "M-phase" from "mixed" of the 1st 
order transition, can thus still be called "the M-phase", 
now from "magnetic". While the QGP at T > 1.4T C is in 
fact dominated by quarks and gluons, it has been shown 
in [28] [31] (in classical and quantum settings respectively) 
that in a mixture of both electric and magnetic plas- 
mas the electric-magnetic scattering dominates transport 
cross sections and viscosity till very high T, thus provid- 
ing the natural explanation for unusually small viscosity 
observed at RHIC and predicting its value at the LHC. 
Its main reason is the so called "caging" phenomenon, 
induced by the Lorentz force. 

The most important consequence of these development 
for the present paper is that in the M-phase of the col- 
lision there is dense magnetic plasma, in which mag- 
netic field is well screened while electric one remains 
unscreened for a significant time. This leads to the cen- 
tral idea of this paper: that QGP produced in central 
part of RHIC collisions should have a "dual corona" in 
which electric (rather than magnetic) fields coexist with 
the plasma, affecting both the overall expansion and the 
propagation of perturbations. We will suggest to use 
"dual magnetohydrodynamics" (DMHD) for the descrip- 
tion of diffuse electric fields in the M-phase, in particular 
study their effect on the velocity of propagation of small 
perturbations. We will further argue below that like solar 
corona, that of QGP should have met ast able flux tubes, 
although microscopically thin ones which cannot be di- 
rectly described by DMHD approximation. 

The word "corona" used here comes from the physics 
of the Sun. Let me briefly remind the reader that it was 
started by Galileo Galilei, who in 1612 spent some time 
observing the motion of the black spots on the Sun and 
correctly concluded from motion of the spots that they 
must resign on a surface of a rotating sphere: he thus ar- 
gued the spots were not shadows of some planets passing 
in front, as it was thought of before. In due time relation 
between the spots and solar magnetism was understood: 
modern telescopes allows one to see the fine structure 
of solar spots, resolving individual magnetic flux tubes. 
Better understanding of solar magnetism came with the 
advance of plasma physics in 1940's and development 
of MHD, which explained both the influence of diffuse 
magnetic field on plasma and formation and mechanical 
stability of the flux tubes. The MHD flux tubes are sup- 
ported by the electron current, while the positive charges 
- the ions - are heavy and dont move. Since it is not a su- 
percurrent, there is inevitable friction and thus metasta- 
bility of the flux tube solutions. 

As at RHIC the central part of the produced fireball 
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FIG. 1: A snapshot of unscreened electric (dual-magnetic) 
field in the M (near-T c ) region of the fireball. Fig. (a) and 
(b) correspond to full RHIC energy and the reduced energy 
(analogous to SPS). 



reaches relatively high temperature T ~ 2T C , we expect 
both E, B fields to be effectively screened there, see the 
central cylindrical part of Fig [I] marked QGP. But in the 
near-T c region vanishing electric component leads to van- 
ishing electric screening mass. This means that plasma 
in the outer cylindrical part of Fig {!] marked M (mixed or 
magnetic) is nearly pure magnetic. It is very important 
to emphasize that although this region on the phase dia- 
gram is represented by a very narrow strip \T — T c \ <T C , 
it corresponds to more than order of variation of the en- 
ergy or entropy density, and the corresponding space- 
time volume in the expansion of the fireball is by no 
means small. I A snapshot of the geometry of the M 
region at some early time is shown in Fig{lj here are 
unscreened electric fields (thin lines) and magnetic flux 
tubes (think lines). The lower plot show similar snapshot 
at collision energy much smaller than at RHIC, planned 
to be investigated in a specialized run. 



B. New structures observed in two and three 
particle correlations 

Three different correlation phenomena have been dis- 
covered in heavy ion collisions at RHIC: 

(i) the so called "cone" [13 [18], is a two-peaked structure 
seen in azimuthal distribution of hadrons on the "away- 
side" from a trigger hadron (the region off quenched com- 
panion jet); 

(ii) the so called "hard ridge" seen on the "same-side" in 
the triggered events [19]; 

(iii) and the "soft ridge" observed in 2-particle correla- 
tions without any hard trigger [20 in the minijet region, 
with transverse momenta p t ~ 1 — 2 GeV . 




FIG. 2: A sketch of the transverse plane of the colliding sys- 
tem: the "spots" of extra density (a) are shown as black 
disks, to be moved by collective radial flow (arrows). Naive 
sound expansion (b) would produce large-size and small am- 
plitude wave: yet the correct solution includes also brighter 
secondary wave (c) of smaller radius. 



(i) The "cone" has been discovered in the 2-particle 
azimuthal correlations like the one shown in Figj3] One 
can see from this figure the disappearance of the "away- 
side" peak at Acj) = tt and appearance of new peaks at 
completely different angle, as one moves from periph- 
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FIG. 3: A set of two-particle correlators from PHENIX col- 
laboration, as a function of azimuthal angle difference A<j>. 
Six pictures are for different centrality classes, indicated by 
percentage of the total cross section. While the most periph- 
eral collisions (the right-lower corner) shows two peaks, one 
near zero and one near <j) — 7r, as in pp collisions, others show 
a minimum for that "away-side" angle and a peak shifted by 
a large angle from it (vertical solid lines). 



eral to central collisions. After discovery of those effects 
there was extensive studies of the 3-particle correlations 
as well. This is a rather complicated subject to go into 
here, let me just say that they has confirmed the observed 
structure is indeed cone-like, and not e.g. a reflected jet. 

(ii) the hard ridge is also seen in 2-particle correlators, 
but plotted on the two-dimensional A<p — At] plane, the 
differences between the azimuthal angles and pseudora- 
pidities of the two particles. The jet remnants make a 
peak near Acj) = 0, A 77 = 0, which was found to sit on 
top of the "ridge" , with comparable width in A<p but very 
wide width A 77. For plots and various features one can 
consult the original talk by Putschke [19 j. Later it was 
shown by PHOBOS collaboration [22] that the rapidity 
range of the ridge extends at least up to \r]\ « 4. 

(iii) the "soft ridge" is found by STAR collaboration 
[20] |2T] without a trigger, in the 2-particle correlations. 

For many experimental details and phenomenological 
considerations related to these phenomena the reader 
may consult e.g. the talks at recent specialized work- 
shop [36] . 

We will return to these observations below, turning 
now to their suggested explanations: 
(i) Stoecker et al, as well as Casalderrey, Teaney and 
myself [30 have proposed that the energy deposited by 
a quenched jet goes into two hydrodynamical excitation 
modes, the sound and the so called diffusion or wake 
modes. The sound from the propagating jet should thus 
create the famous Mach cone, in qualitative agreement 
with the conical structure observed. 



(ii) One early model for "hard ridge" has been introduced 
in my paper [24] . It relates it with the forward-backward 
jets accompanying any hard scattering, providing extra 
particles ( "hot spot" ) widely distributed in rapidity. This 
idea is then combined with the one suggested previously 
by Voloshin [23] . namely that extra particles deposited 
in the fireball would be moved transversely by the ra- 
dial hydrodynamical flow, should produce a peak at cer- 
tain azimuthal angle corresponding to the position of the 
hot spot, see Fig(2|a). While particles of the ridge are 
separated by large rapidity gaps and cannot communi- 
cate during the expansion process, their azimuthal emis- 
sion angles remain correlated with each other because 
they originate from the same "hot spot" in the trans- 
verse plane. 

(iii) Similarly, transverse hydro boost of "hot spots" was 
used for the explanation of the "soft ridge" by McLerran 
and collaborators [25] [26] . They have pointed out that 
the initial state color fluctuations in the colliding nuclei 
would create longitudinal "color flux tubes" without any 
hard collisions. As these tubes are being stretched be- 
tween two fragmentation regions of the colliding nuclei, 
they also lead to long-range rapidity correlations. 

So, at a very qualitative level the origin of all three 
phenomena seem to be explained: yet at more qualita- 
tive level a lot of puzzles appear. As an example, con- 
sider the simplest of them, the "soft ridge" . As discussed 
in [25] [26] . the initial stage (proper time r ~ l/Q s ~ 
0.2 fm/c where Q s ~ lGeV is the so called saturation 
scale at RHIC) can be discussed using classical Yang- 
Mills equations: thus color fluctuations naturally appear. 
However, the observed pions come from final freezeout 
time, separated from the initial "glasma" era by much 
longer time r ~ 10 fm. This is certainly so, as the ex- 
planation heavily relies on radial hydro velocity and thus 
it has to wait till the hydro velocity is being created. As 
we will argue below, there are many reasons why one 
might have expected nearly complete disappearance of 
this signal during this time. 

(We should note that one recent hydrodynamical study 
of the fate of tube-like fluctuations, by Takahashi et 
al [27], report observing both the ridge and the cone- 
like structures from initial fluctuations without any jets. 
There is however no explanation in the paper of why this 
may be happening, so we have to wait for further works.) 

Common to all three cases is deposition of some addi- 
tional energy (or entropy), on top of the "ambient mat- 
ter". The number of correlated particles in all of them 
constitute a small (~ 10 -3 ) fraction of the total mul- 
tiplicity: thus they can only be seen in a high-statistics 
correlation analysis. Furthermore, simple estimates show 
that it would not be possible to detect any trace of that 
tiny perturbation if it would be distributed over a signif- 
icant fraction of the fireball: the only possibility is that 
it remains well localized in transverse direction. 

Smallness of perturbation in respect to total system 
size by itself does not guarantee that the perturbations 
is small locally, in respect to local density of ambient 



matter. However it will become so if perturbation would 
give rise to divergent conical (or cylindrical, or spheri- 
cal) waves, see Figj2|b). Similar to circles from a stone 
thrown into a pond, initial perturbation may become 
some waves, with basically nothing left at the original 
location at later time. Even without dissipation, ideal hy- 
drodynamics predicts that the final radius of those waves 
is given by the "sound horizon" 
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As we will detail below, by the the freezeout proper time 
Tf ~ 10 — 15/ra/c, this distance is not small, ~ 6/ra 
or so, since the speed of sound changes between c s = 
l/\/3 ~ .58 in QGP and about .3 at its minimum near T c 
The amplitude of the wave is decreasing accordingly, and 
the width of <j) distribution grows, making us wandering 
if any trace of the perturbation can remain observable. 

And yet, we do observe all three correlations, as if 
nothing happened to them during rather long time of 
the hydro process, ~ 10 fm/c. This is the puzzle dis- 
cussion/resolution of which is the main objective of this 
paper. The idea behind it is that there can be some 
reasons providing second unusual mode of propagation, 
with reduction or maybe even vanishing of the speed of 
its spread and lead to an observable structures at late 
time, see see Figj2^c). 

(In the case of a cone, additional consideration is that 
the "wake" mode, behind the jet, is - in contrast - not 
expanding or weakening: and yet it is not observed. In 
the case of ridges, large size of waves comparable to nu- 
clear radius will make the radial flow directions be rather 
different at different places, widening the peak in azimuth 
well beyond what is actually observed. The puzzle is es- 
pecially clear in the case of "ridges" , whose explanation 
heavily rely on substantial hydrodynamical flow velocity, 
which cannot be formed promptly and is known to be 
developed only by the freezeout time.) 

Let me now add few more important observations on 
the soft ridges. The spectra of particles in cone and 
ridges, as well as their composition (not shown) are 
drastically different from jet remnants [19] . Particularly 
telling is large baryon/meson ratio which clearly indi- 
cate that their existence is related to the ambient mat- 
ter boosted by the hydrodynamical flow. The boosted 
baryons and sharpening of the <p peak nicely confirm that 
the particles of the ridge do come late, from the final 
freezeout of the system. 

Further confirmation of hydro origin of ridges comes 
from the centrality dependence of the angular width of 
the ridge: the peak in azimuth sharpens for more cen- 
tral collisions, see Figj4j This happens because of two 
interrelated effects, both well documented. For central 
collisions there is (i) an increase of the radial hydro ve- 
locity, accompanied by (ii) a substantial decrease in the 
freezeout temperature (which goes from T c « 170 MeV 
in peripheral down to Tf w 90MeV for central collisions. 
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FIG. 4: Azimuthal width of the ridge as a function of central- 
ity, from the talk of L.Ray [36] . The red squares and black 
dots are for 200 and 62 GeV Au Au collisions. Green trian- 
gles correspond to the minijet model without matter effects 
shown for comparison. 



The absolute magnitude of the peak's width also pro- 
vides a limit on how widely the "spot" can evolve, from 
its creation at inital time to the freezeout moment. In 
Figj5]we have plotted the shape of azimuthal peak pro- 
duced by (semi) circles of radii 1..6 fm. We have super- 
imposed Hubble flow and calculated this angular distri- 
bution. As one can see from this figure, the width of 
the distribution grows - it is 0.57, 0.56, 0.69, 0.76, 0.83, 
0.89 for the 6 curves shown. Moreover, the distribution 
shapes become very different from that observed, with 
two maxima shifted from <p = (corresponding to di- 
rection of flow at two points at which the circle intersect 
the fireball boundary). Comparing such distribution with 
observations, e.g. their width with those shown in Figj4] 
we concluded that for central collisions the radius of a 
spot at freezeout is restricted to be R(rf reezeout ) < 3 fm 
or so. As we already argued in the Introduction, this is 
by about factor two smaller than the radius of the sound 
horizon with realistic speed of sound . 

Having mentioned the agreement with overall hydro 
- radial and elliptic flows - it is only fair to mention 
their important differences. The (soft) ridge has dra- 
matic centrality dependence shown in Fig j6j (The lines 
and shaded area marks GLS is some simple scaling ex- 
pected from noninteracting minijet event generator: it 
only describe the data at peripheral bins at small den- 
sities). Furthermore, note that the comparison of the 
200 and 62 GeV AuAu data shows that the transition 
point seem to be at the same transverse particle density 
Pc — (3/2)dN c h/dy / S « 3. One might naively think that 
for p < p c the matter is simply too dilute to show hy- 
drodynamical effects: but in fact both radial and elliptic 
flows (e.g. V2A2) are quite smooth at such centrality and 
show no rapid changes. 
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FIG. 5: The dependence of the visible distribution in az- 
imuthal angle on the width of the (semi) circle at the time 
of freezeout. Six curves, from the most narrow to wider ones, 
correspond to the radius of the circle 1,2,3,4,5,6 fm, respec- 
tively. The original spot position is selected to be at the edge 
of the nuclei. The distribution is calculated for particle of 
p t = 1 GeV and fix freezeout T f = 165 MeV. 



This difference between ridges and overall hydro flows 
may be directly related to the main dilemma of this pa- 
per. If the cones and ridges are hydrodynamical, then 
why can they be so different from overall hydrodynami- 
cal expansion in their centrality dependence? 
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FIG. 6: The dependence of the amplitude and azimuthal 
width of the "soft ridge" on the transverses multiplicity den- 
sity, from the talk of L.Ray [36]. Closed and open points are 
for 200 and 62 GeV AuAu collisions. 



Strong temperature dependence can in principle be 
related to the issue of timing of the energy deposition. 
There is a difference between timing of the "cone" and 
"ridges": while the latter obviously originate early, for 



"cones" the exact energy deposition time/place depends 
on the jet quenching mechanism. As gluon (or light 
quark) jets move with a speed of light, by the time of 
the order of nuclear size ~ 6/m/c they either leave 
the fireball, or are already completely quenched. Re- 
cently purely geometrical study of angular distribution 
of quenching [34] indicated that most of jet quenching 
rate should happen in the near-T c region. It is surpris- 
ing, taking into account much higher density of QGP at 
earlier time. 

Another dramatic finding, pointing to the same direc- 
tion, was unexpected observation of similar "cones" by 
CERES and NA49 collaborations at much lower collision 
energy of CERN SPS (see recent summary by Appel- 
hauser in [36). Since at the SPS energy there is practi- 
cally no QGP phase, it can only be there starting in the 
near-T c (mixed) phase. 

Obviously we would like to see what happens with the 
ridges at lower collision energies. Note that both ex- 
planations we propose in this work have problems with 
QGP away from T c : it is hard to stabilize the flux tube 
there and also impossible to stop expansion of the sound 
waves with rather large sound speed c s = 1/ y/3. And yet 
ridges disappear in very peripheral collisions: we would 
like to know what happens as the collision energy gets 
lower. (Those questions are presumably be addressed by 
the expected scan down in RHIC energy, planned in the 
nearest future.) 

Another surprising experimental fact is quite large 
value of the cone angle, deduced from 2 and 3-particle 
correlators. It seems to be in the range 9m = 1.2 — 1.4 ra- 
dians (not too far from ir/2 = 90° or cylindrical waves!). 
The Mach formula gives the speed of pertinent perturba- 
tion to be about 

< Vwave >= C0S(6 M ) ~ 0.2 (1.3) 

well below the expected speed of sound (except maybe 
near T c ). So again, it is either (A) a coil effect, reducing 
expansion, or (B) a sound with a nontrivial production 
mechanism. 



C. Introductory information about the fireball 
hydrodynamics 

As we already mentioned in the Introduction, hydro- 
dynamical description of the explosion following heavy 
ion collisions - known also as a "Little Bang" - are by 
now widely accepted as a good quantitative description 
of the process. There is no point to repeat here how they 
describe radial and elliptic flow of matter. However we 
think we need to remind the reader the geometry and 
timing of this process at RHIC energies, with proposed 
simplification of the flow, to be used below. 

Although the "Little Bang" is short in absolute time, 
it is a complicated process in which the system undergoes 
expansion, cools down and go through different QCD 
"eras" and finally freezeout into the observed hadrons 
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The value we will use in the examples shown below is 
H « 0.08 /m _1 , as appropriate for RHIC hydro. 

(Let us also remind in passing that a lot of misun- 
derstanding in the years of 2001-2004 of hydrodynamical 
results were related with incorrect determination of the 
freezeout surface. By no means it corresponds to the the 
fixed-proper-time r = const surface, or an isotherm (e.g. 
T = 120 MeV shown in Figj7|, as several hydro groups 
incorrectly assumed. That wrong assumption had lead to 
several dependences on collision energy, rapidity and cen- 
trality of the elliptic flow in disagreement with data. The 
freezeout is equality of the reaction rate and the expan- 
sion rate: and the latter depends on absolute system size 
and particular point: in fact it is closer to the upper-left- 
to-lower- right diagonal of this figure. If this is properly 
taken into account - e.g. by switching to hadronic cas- 
cades as in pQ - all these problem go away, for details see 

e.g. m) 



FIG. 7: The hydro-dynamic solution for RHIC collision en- 
ergy, the horizontal axes is the transverse radius r (each large 
tick 1 fm), the vertical is proper time (each small tick 1 fm). 



II. ACOUSTICAL WAVES IN EXPANDING 
FIREBALL WITH VARIABLE SPEED OF SOUND 



at proper time r ~ 10 — 15/ra/c. In Fig 7] from [1 thick 
lines show contours of constant energy density. ei2o de- 
notes the energy density where T = 120 MeV. e# and cq 
denote the energy density where the matter shifts from 
hadronic to mixed and mixed to a QGP, respectively. 
(Ut) denotes the mean transverse rapidity weighted with 
the total entropy flowing through the energy density con- 
tours. Walking along these contours, the line is broken 
into segments by dashed and then solid lines. 20% of 
the total entropy passing through the entire arc passes 
through each segment. The thin lines show contours of 
constant transverse fluid rapidity (vt = tanh(?/ T )) with 
values 0.1,0.2,. ..,0.7. 

The main lessons from this plot are: 

(i) The times spent in three phases (QGP, mixed M 
or near-T c , and hadronic H) are all comparable at the 
center of the fireball, r = 0, while the QGP phase 
disappears toward the rim of the Au nuclei r > 6.5 fm. 

(ii) In the QGP era the flow is building and is therefore 
complicated. However in the later mixed and hadronic 
eras the flow picture substantially simplifies, as is seen 
from the fact that the thin lines (of fixed transverse ra- 
pidity) are nearly exactly vertical and nearly equidistant. 

This second feature we will use to simplify the picture, 
using by the so called Hubble flow 



v(t, r) = Hr 



(1.4) 



with the Hubble constant H = const(r,t) being nearly 
space and time-independent. If so, expansion is quite 
simple 



A. Model equation of motion 

Now we turn to discussion of the evolution of small 
perturbations sitting on top of overall (Hubble-like) ex- 
pansion. Equations for this case have been worked out by 
Casalderrey-Solana and myself in Ref. [8 , and applied 
to "conical flow" [30 from quenched jets. In this paper 
we have already found that in the vicinity of the QCD 
phase transition there is a wave splitting phenomenon. 
The framework we will use to study the effects of the vari- 
able speed of sound and matter expansion we have looked 
for the simplest example possible, keeping the problem 
time-dependent but homogeneous in space. This can be 
achieved in a Big-Bang-like setting in which the space is 
created dynamically by gravity. Consider a liquid in flat 
Freedman-Robertson- Walker metric : 

dr 2 = dt 2 - R(t) 2 [dr 2 + r 2 (d0 2 + sin 2 6d(j))] (2.1) 

where the parameter R(t) (the instantaneous Hubble ra- 
dius of our "universe") is treated as external (not to be 
derived from Einstein equations but from hydrodynam- 
ical solution such as shown in Figj7]). For isotropic ex- 
pansion only the longitudinal projection u^Tf^ is needed, 
which is the equation of entropy conservation leading to 



dt 



(s(t)R 3 (t)) = s(t)R(t) 3 = S 



(2.2) 



r(t) = r(0)exp(Ht) 



(1.5) 



for any R(t), provided expansion is adiabatic. 

A simple substitution of a variable c s (t) into the equa- 
tions of motion for perturbations is inconsistent. One 
should instead find a correct non static solution of the 
hydrodynamical equations and only then, using this solu- 
tion as zeroth order, study first order perturbations such 
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as sound propagation. The linearized equations for hy- 
drodynamical perturbations in this background has been 
derived in [8]. Using the normalized perturbation 



Profiles starting at r=0 during the QGP phase at tau=0..5 
0.4- 



R A 5T 



-oo 



(2.3) 



one may eliminate other components of the stress tensor 
and get the following single equation 



(2.4) 



f dt [(Sc 2 s (t) - l)d t R(t)] + (3c 2 s (t))^^(d t e) = 0. 



R(t) 



R(t) 



In the derivation we have not assumed any particular 
expansion function R(t) or particular equation of state, 
just general thermodynamic relations. 

In Ref. [8 we have used a bit different time vari- 
able and Fourier decomposition in space, reducing the 
problem to an oscillator with the time-dependent fre- 
quency and specific exciting force (the third term) which 
is negative for c 2 s < 1/3. Note that it creates amplifica- 
tion of dimensionless perturbation, which is similar to the 
case of Universe expansion effect on small perturbations, 
running away from each other. 



B. Generation of the secondary wave 

Let us remind the setting in which the solutions were 
studied. We have already explained that we expect the 
expansion to be exponential R(t) = exp(Ht) and the H 
value was fixed from transverse matter expansion. For 
simplicity, we use the same expansion in 3d, although in 
heavy ion collisions the longitudinal expansion is differ- 
ent. 

The main ingredient is the variable speed of sound 
c s (t). As we have already emphasized in the Introduc- 
tion, at early stages at RHIC the matter is believed to be 
in the form of quark-gluon plasma (QGP), and thus wit h 
Cq GP ~ 1/3. This makes the third term in our eqn (2.4): 
which is nice since in this stage the Hubble flow is not yet 
a good approximation. In the near-T c region the energy 
density is increasing much more rapidly than the pres- 
sure, which makes matter "soft" and c 
minimum, known as the "softest point 
is rising again, in the hadronic "resonance gas". Our 
model-dependent time histories at two transverse posi- 
tion r = 0, 6/m are depicted in Figj9j for three different 
scenarios we would study. 

Solutions for eqn ( |2.4[ ) corresponding to Figj9^a) (the 
sharpest change) for initial Gaussian perturbation is 
shown in the Fig(8j In each cases we show in three pic- 
tures subsequent stages of evolution for r * e(t, r) as pro- 
files taken every fm/c. We start with 5 curves of the 
fireball history, Figsj8ja) corresponding to the QG P era. 
Since c 2 s « 1/3 the main new (third) term in (2.4) is ab- 
sent, and our solution corresponds to "naive evolution" of 
the original perturbation into expanding primary wave. 



2 dropping to its 
Aftre that c s (t) 




Profiles starting at r=0 during in the mixed phase at tau=6. 
.10 




Profile starting at r=0 during the hadronic phase at tau=ll. 
.15 




FIG. 8: The snapshots of the evolution of perturbation in 
the QGP (a), mixed (b) and hadronic (c) phases 



No visible trace of the perturbation remains at the initial 
position, and their are no secondary waves. 

In Figsj8^b) we show what happens during the mixed 
phase: the evolution of perturbations slows down as ex- 
pected. One also finds growth in amplitude, and also the 
secondary wave starts to appear. Figs|8jc) show that in 
the hadronic phase both waves start propagating outward 
- if there is time left to freezeout. 
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FIG. 9: (Color online) The sound velocity c s as a function of 
proper time r (in fm/c), for the fireball center (r=0, red) and 
the rim (r=6 fm, green). (a,b,c) are three variants of the time 
dependence, corresponding to the results display in the next 
figure. 



Final profiles at r=0 (black) and r=6 (brown) 




Final profiles at r=0 (black) and r=6 (brown) 




FIG. 10: A comparison between the final profiles, in the two 
cases studied. The freezeout is assumed to be at 15 fm/c 
at the center, r=0, and 11 fm/c at r=6. (a,b,c) are three 
variants of the time dependence shown in the precious figure. 
Note how the amplitudes of the secondary waves gets smaller 
as the time variation of speed of sound gets more smooth. 



We have further studied how the amplitude of the sec- 
ondary wave depends on the variation of speed of sound. 
For this purpose we have used e.g. 3 scenarios depicted 
in Fig. [9j In Fig 10 (a-c) we show the final profile of 
the wave at freezeout, for two histories corresponding to 
r = (the fireball center) and r = 6 fm (fireball rim). 

Let us start with (a) in which one finds re of the sec- 
ondary wave at the rim of about the same magnitude 
as the primary wave. It means that the sound intensity 
(which scales as 1/r 2 ) is for the secondary wave about 
an order of magnitude larger. This would make the sec- 
ondary wave much more likely candidate to be observ- 



able. 

Furthermore, the azimuthal angular width of the ridge 
observed experimentally is less than 1 rad, see Fi^4] 
Given the magnitude of the rapidity of the radial flow 
- up to 0.7 near the edge of the fireball, as seen in next- 
to-last thin line in Fig. [7]- and known thermal spread 
of secondaries at freezeout, one can explain the observed 
width assuming the spot size remains small as compared 
to the fireball radius. This condition is fulfilled only for 
the secondary wave, while for the primary one one should 
include extra spread of the radial flow directions. Thus 
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the secondary wave offers possible resolution of the puz- 
zle. 



As one can see from the solution shown in Fig. 10 
(b,c), the effect is unfortunately not quite robust: the 
amplitude of the secondary wave decreases if the speed 
of sound is changing more smoothly. 



III. DUAL MAGNETOHYDRODYNAMICS 
(DMHD) 



Thinking about possible applications of DMHD, one 
should first clearly separate two opposite limits: (i) the 
weak field case, withp^ e ^ <C Ppiasma] an d (h) the strong 
field case, with p field ^ Ppiasma- We will discuss them 
subsequently. In the former case matter properties is 
only weakly affected by imbedded field, while strong field 
region expands till the pressure balance is reached, ex- 
pelling plasma from it. The situations with weak "dif- 
fuse" fields as well as strong fields, creating flux tubes 
with no plasma inside, is well known in e.g. solar plasma. 



The major new idea we put forward in this work is 
the suggestion to use the dual Magnetohydrodynamics 
(DMHD) as an approximation effectively valid in the 
near-T c (or M) region. This proposal would of course 
lead to multiple consequences, of which we will discuss 
only the simplest ones. Qualitatively, the main effect of 
the electric (dual-magnetic) fields in plasma can be in- 
corporated simply by including the fields stress tensor 
together with that of the plasma. It would lead to ex- 
tra "elasticity" (pressure), helping the overall expansion 
a bit. Furthermore, as fields are directed longitudinally, 
one gets certain anisotropy of the medium, with the speed 
of perturbation depending on its angle relative to the 
beam (and field) axis. 

The notations we are going to use are simply dual to 
standard ones in MHD. Thus the "magnetic current" of 

moving monopoles would be denoted by j, the coupling 
constant g, and the (gluo)electric field E related to dual 

H, but with different normalization. The pressure of the 
field is 



E 2 



IP 

8tt 



(3.1) 



where we keep different 4tt in normalization, Thus, apart 
of tildas, Maxwell eqns look familiar, same as in text- 
books. In the infinite conductivity limit, those can be 
written as two dual-magnetic eqns 



divH 
= curl[vH] 

complemented by Euler eqn of hydrodynamics 



(3.2) 



OH 

~dt 



dpvi 
~dt 



dx k 



(3.3) 



1 



U lk = P ViV k +p5 ik - j- (HiHk - (l/2)H 2 S ik ) (3.4) 

appended by the magnetic stress tensor, as well as the 
usual matter continuity eqn 



dp dpvj 
dt dxj 



= 



(3.5) 



A. Perturbations in the (dual)-magnetized plasma 

In this work we restrict ourselves to the simplest prob- 
lem, that of propagation of small- amplitude waves, in the 
simplest geometry, in which only one (longitudinal or z) 
component of the field is nonzero. Assuming the field 
to be permanent and homogeneous of some amplitude 
Bo = const(x, t) one can linearize the MHD equations, 
look for plain wave solutions and determine the disper- 
sion relation for the waves. 

Qualitatively it is not hard to tell what is going to 
happen: depending on the relative sign of the pressure 
perturbation and that of the field, there would be two 
solutions, one with a speed large and one smaller than the 
sound speed. This is well documented textbook problem, 
see e.g. chapter 69 of [7j, so we just mention the final 
expression for the velocities of two "magnetosounds" or 
Alfven waves is 



#2 
47TC 



(3.6) 



47T6 



2\2 



H 2 cos 2 0c 2 



-,1/2 



where c s , e is the speed of sound in "unmagnetized" 
medium without field and the energy density, 9 is the 
angle between the field strength and the the direction of 
the wave propagation. Note a case in which the wave 
goes transverse to the field ( cosO = ) in which the 
lower mode has zero speed. 

In the case of jet quenching - when the jet direction is 
more or less up to experimentalist to pick - general shape 
of these waves can be complicated. However when the jet 
(or original charge fluctuation) propagates longitudinally, 
in the same direction as the field, the problem is axially 
symmetric and results in general in two cones. The angles 
of their propagation can be obtained from the previous 
expression, in which the l.h.s. is substituted by Mach 
relation u — > cos(6)v, v is the velocity of the jet, and 
solve it for the cosO. The resulting equation can be solved 
analytically, giving 



cosO ■ 



Ml 

47T6 



(v 2 



(3.7) 
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FIG. 11: (Color online) The solid (red) lines in figs (a) 
and (b) are for the velocity of the outer and inner cones, 
respectively. The dashed (blue) lines are v 2 * z 2 versus 
z = cos(0). The crossing of lines - correspondence of the 
Mach condition to the wave velocity - give the angle of the 
conical waves. All curves are for c 2 = 0.09, which is the 
sound velocity at T c according to the lattice data. In (a) 4 
solid curves are for H 2 /4tt€ = 0.1,0.2,0.3,0.4, top to bot- 
tom, and 4 dashed for v 2 = 0.1,0.2,0.3,1., from bottom 
In (b) H 2 /4tv€ = 0.1,0.2,0.3,0.5 for solid lines and 



up 

„,2 



0.05,0.06,0.07,0.08,0.09. 



or zero cosO = or no solution. One obvious condition 
is that a jet should be "supersonic" v 2 > 

Better insight into this equation is provided by Fig|TT| 
which display several values of the field energy density 
relative to matter and jet velocity v for which nontrivial 
cone angles exist. As one can see from the lower plot, 
only for rather slow jet there are solutions with nonzero 
cosO: otherwise the cone solution is at cosO = corre- 
sponding to zero velocity or non-expanding (stabilized) 
field region. In general, note that the scales on two figures 
are different: the angle of the outer cone is larger than for 
ordinary Mach cone, while the inner one is small if not 
zero. Therefore, qualitatively we return to the picture 
depicted in Figj2^c). For obvious reason, one may think 
that the slower waves are brighter: perhaps we observe 
those. 



B. Flux tubes 

In the previous subsection we considered the electric 
field to be constant in space. If instead we have a spot of 
such field, localized in transverse plane, inevitable there 
is nonzero dH z /dr, which is a part of curl(H) and by 
Maxwell equation it should be proportional to (dual) cur- 
rent j0. This tells us that a flux tube solution must have 
a "coil" with a current running around and trying to can- 
cel the field outside the spot. Ideal DMHD has axially 
symmetric solution - the classical flux tube - similar to 
what is used in solar physics. Two nonzero equations for 
r-dependent p(r)j^(r),H z (r) of the set are 



dp(r) 
dr 



1 dH z (r) 
47r dr 



(3.8) 



(3.9) 



from which it follows that the pressure is balanced in a 
simple way 



p(r) 



H 2 z (r) 
8tt 



= p(r = oo ) 



(3.10) 



Two equations for three functions mean that there is 
functional freedom to select the tube profile. Dissipa- 
tive terms proportional to viscosity and (1 /conductivity) 
can also be accounted for: they lead to more equations. 

Furthremore, MHD flux tubes in electromagnetic plas- 
mas is different from that of the plasma of monopoles 
under considerations: unlike electrons and ions, the 
monopoles and antimonopoles of opposite charges have 
the same mass. Thus the monopole current actually con- 
sists of two components of different charge, counterrotat- 
ing in the opposite directions. If the magnetic plasma 
is sufficiently strongly coupled, their mutual rescatter- 
ing would produce significant friction and short lifetime 
of such configuration. However the main distinction be- 
tween the macroscopically large flux tubes, which can be 
described by MHD equations, and the QCD flux tubes is 
that the latter have microscopic transverse size, smaller 
or comparable to the mean free path in the medium. 

Let us at this point remind the reader brief history of 
flux tubes in QCD, also known as QCD strings. Early 
ideas that such flux tubes are surrounded by a magnetic 
supercurrent due to presumed "dual superconductivity" 
in the QCD vacuum [39] were refined in magnetic effec- 
tive theory [9] and confirmed by multiple lattice studies 
such as [10]. With the advent of "magnetic scenario" 
[31] [32] for the near-T c region, it was suggested that "nor- 
mal" (Bose-uncondenced) magnetic quasiparticles would 
also be able to create a "coil" around electric flux tubes, 
sufficient to stabilize them, above certain density. 

Stability condition of microscopically small metastable 
flux tubes, created by scattering of monopoles on the 
electric flux, has been worked out in two papers by 
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Liao and myself [35] [38]. If the flux tube is small and 
monopoles can penetrate into it, their contribution to the 
positive or negative current depends on the (cylindrical) 
partial wave: as a result expression for the current j are 
rather involved. There is no need to describe these cal- 
culations here: let me just quote the final condition for 
the mechanical stability of the flux tube 



7rM 2 c 2 k B T 
h 2 g 2 (T)n{T) 



1/2 



< 0.13 



(3.11) 



where M(T),T, g(T),n(T) are the monopole mass, the 
temperature, the magnetic coupling and the monopole 
density, respectively. The numerical value in the r.h.s. 
follows from numerical solution for the flux tubes subject 
to quantum scattering by magnetic monopoles. 

So, what is the T range in which the density of both 
condensed and "normal" monopoles is sufficient to sup- 
port the flux tubes? Are there any phenomenological or 
numerical evidences that such flux tubes actually exist? 
In Refs[35 , 38 the main input idea was based on lattice 
data of the interquark potentials at finite T. In brief, the 
central observation is large difference between the free 
energy F(T, r) and the potential energy 



V(T,r)=F(T,r)+TS(T,r) 



(3.12) 



associated with quark pair at distance r. We will not 
show the potentials themselves but just their effective 
string tensions for both, calculated as a slope of the linear 
part extracted from [43] , is shown in Figfl2| The phys- 
ical difference between the two, first discussed by Zahed 
and myself [40] in the context of hadronic spectroscopy 
at finite T, is that F corresponds to adiabatically slow 
motion of the quarks, slow enough to produce maximal 
entropy possible and reach thermal equilibrium at any r. 
However when quarks are moving with certain velocity 
v away from each other, only a fraction x of the maxi- 
mal entropy can be produced (because of Landau- Zener 
argument on level crossing): thus the effective potential 
would be 

V eff (T, r) = F(T, r) + (1 - x)TS(T, r) (3.13) 

For relatively rapid motion in which no entropy is pro- 
duced, x = 0, and one returns to V(T, r). These argu- 
ments have been important to discussion of charmonium 
survival at RHIC as well as say dominance of baryons 
[33] in the near-T c region seen on the lattice. 

The difference between the two potentials is discussed 
in detail in [38]: in short F(T, r) has been related to "su- 
percurrent coil" and V(T, r) to "normal met ast able coil". 
If so, large peak at T c of the tension of V(T, r) is related 
[35] [38] to a peak in "normal" monopole density at T c . 
The condenced and normal monopole density needed to 
accomodate flux tubes with such tensions have been de- 
rived in Refs[35 , 38 and compared with direct lattice ob- 
servations of monopoles. The summary of those studies 
is that metastable flux tubes seem to exist at T < 1.4T C , 
changing from stable to metastable around T c . 
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FIG. 12: Effective string tensions in the free energy <jf(T) 
(squares, from [43 ) and the potential energy ov(T) (dia- 
monds, extracted from 1411) 



C. Production and breaking of the flux tubes 

Since the cones and ridges are all created by electri- 
cally charged particles propagating in the plasma, they 
do start with rotating magnetic field related to "dual 
Faraday" effect shown in Fig|l3] While in electric plasma 
the acceleration happens radially, the magnetic objects 
are accelerated by the magnetic field. For a relativist ic 
charge those are circles concentrated in a pancake of the 
width St ~ I/7 ~ y/\ — v 2 with the spatial distribution 



B 



jg e [vf ± ] 



[j 2 (z-vt) 2 + r 2 ± ] 3 / 2 



(3.14) 



Monopoles (with magnetic charge g m ) experience instant 
kick, in the corresponding direction. Its magnitude has 
no gamma factor and is proportional to the product of 
electric and magnetic couplings g e g m /4:7r = n, a Dirac in- 
teger. (For elementary quark/gluon electric charge and 
for Polyakov-t'Hooft monopoles it is just ). Of course 
electric charge fluctuations can lead to larger charge val- 
ues as well. This velocity (near-) instantaneously pro- 
duces a current (or a "coil') running around the flux. 
The questions under which condition it is robust enough 
to contain the field into a static flux tube will be discussed 
elsewhere. 

The interpretation of the linear part of the potential 
energy coming from the lattice V(T, r) as a metastable 
flux tube mentioned in the previous subsection, Fig|l2] 
can be used to provide estimates of the absolute amount 
of energy/matter involved. At T c the values are astonish- 
ingly large: the energy per unit length (tension) reaches 
at its peak 5 GeV/fm, corresponding to the entropy den- 
sity of about 30/m _1 . One may ask if those parameters 



13 



can be compatible with observations of the ridges/cone 
at RHIC. 

Let me start addressing this issue from the point of 
view of energy first. The total energy of the flux tube, 
or the work done by its tension on the departing charges 
(large- x valence quarks) is SE = ra. For the M-phase 
lasting r ~ 5fm/c and the proposed V-potential tension, 
one gets the energy loss of about 25 GeV. Keeping in 
mind that the total energy of each nucleon in the center 
of mass is 100 GeV, at RHIC energy used for heavy ions, 
and that each nucleon has three valence quarks plus sea 
plus gluons, we conclude that it is comparable to quark 
total energy. So, if the string would not break, it would 
be able to transfer valence quarks from the fragmentation 
region to midrapidity: and we know from experiment 
that it happens with very small probablity. 

The conclusion is then that such strings must get bro- 
ken at time shorter than r ~ 5/m/c used in this esti- 
mate. This is also known to be true for the usual (vac- 
uum or T = 0) strings: for those the rate of breaking 
has been phenomenologically extracted from hadron de- 
cays, especially of hadrons with large angular momenta 
(Regge trajectories) and event generators based on the 
Lund model. 

Naively, one might think that "metastable" flux tubes 
in the M-phase should have smaller lifetime than the 
vacuum QCD strings. First, their decay can be related 
not only to (i) the quark pair production, as in vacuum, 
but also to (ii) just picking up quarks from the ambient 
matter. However, looking at numbers more closely one 
finds that it is not nacessarely so. 

Quark pair production is described by Schwinger 
fermion pair production rate (the leading exponent only) 
in a constant electric field E 



1 which explains small density of these quasiparticles and 
thus small probablity of a string breaking by picking up 
free charges. 



dW 



a 2 E 2 



-exp 



ttMI 
' eE 



(3.15) 



where M is the charged particle mass. The tension scales 
as a ~ E 2 r\, where E, r±_ are electric field and transverse 
size of the flux tube. There is a universal flux, so Er\ = 
const, and thus in any change of the tube the field scales 
as the first power of the tension. If tension grows by 
some factor, e.g. <j(T c )/<j(0) ~ 5 as suggested by lattice 
potential V, the field should grow by the same factor. 
Naively, it would greatly reduce its lifetime. 

Yet the effective mass of quark quasiparticles near T c 
is also quite different from "constituent quark mass" in 
the vaccum (T = 0) M q w 330 MeV . For example, di- 
rect lattice studies [42] show that M qj M g ~ 800 MeV 
at T = 1.5Tc, with perhaps even larger value at T c . 
Thus the combination which enters the exponent of the 
Schwinger formula M 2 /cr is not decreasing but rather 
grows. It suggests a rather counterintuitive conclusion: 
flux tubes at T c may have smaller breaking rate than in 
the vacuum. 

The same argument, based on large quark/gluon mass 
near T c , leads to small Boltzmann weight exp(—M/T) <C 
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FIG. 13: Schematic demonstration of magnetic solenoidal 
by Dual Faraday's law, see text. 



Remarkably, recent RHIC data provided direct experi- 
mental indications for enhanced stabilty of the flux tube 
in matter relative to pp. We will use those from PHO- 
BOS collaboration, which has large rapidity coverage of 
their silicon detector. Fig 14 'a) shows that the number 



of charged particles in a cluster observed in AuAu col- 
lisions (points) is about twice that seen in correlation 
studies of the pp collisions (shaded horizontal region), so 
that they reach the size of 6 charged particles (or 9 to- 
tal). (The exception are central collisions on the right, 
in which case extensive hadronic after-burning kills the 
correlations.) The figure (b) shows that the produced 
clusters are not near-isotropically decaying resonances as 
in pp, (shaded horizontal region), but are instead more 
extended in rapidity. This last fact shows their direct 
relation to the "soft ridge" and flux tubes. Similar CuCu 
data (not shown) demostrate basically the same clus- 
ters, provided the same centrality is taken: this shows 
that geometry and surface-to- volume ratio is important. 
Taken together, we interpret those clustering data as di- 
rect proof of significant changes in the flux tube decay 
parameters in AuAu relative to pp: the tubes apparently 
gets denser and decay less frequently, into larger pieces. 

Let me also try to connect the multiplicity in the clus- 
ters observed with the entropy as seen in lattice poten- 
tials. At the peak the entropy per length is as high as 
dS/dL = 5GeV/ fm/T c ~ 30. If this entropy all goes 
into final pions at freezeout, one can estimate the cluster 
absolute length. Neglecting the pion mass - that is con- 
sidering matter to be ideal massless bose gas at freezeout 
- one find that entropy is 3.6 per particle, or S ~ 36 for 10 
pions in a cluster. Thus if the decay takes place in the M 
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FIG. 14: (Color online) The cluster size (a) and width (b) 
as a function of centrality (cross section fraction) in AuAu 
collisions at RHIC full energy, from PHOBOS [37 . The size 
is the number of charged particles associated with the cluster 
and the width is in rapidity. 



phase, the cluster corresponds to a string length of 1.2 fm 
or so. However if the decay rate is small and the string 
survivies through the M-phase to hadronic phase, the 
string tension is reduced to the vacuum value 1 GeV/fm, 
the length corresponding to the observed cluster would 
be as long as 6 fm. In principle this information can be 
directly related to the r.m.s. rapidity width of the clus- 
ter seen from fig, 14 'b). We are planning to compare it 



with some PYTHIA-like simulations to tell whether those 
clusters can be consistently reproduced. 



D. Discussion and predictions 

In this section we have suggested two consequences of 
the small electric screening mass in the M-phase (near 
T c region): (i) unscreened bulk electric fields and (ii) 
metastable (and possibly even relatively long lived) elec- 
tric flux tube. Taken together, they were called "QGP 
corona", with reference to similar phenomena in solar 



plasma. We had tried to connect this picture with phe- 
nomenology of fireball perturbations, namely ridges and 
cone. Now we would like to speculate further along 
this line, pointing out some further consequences of the 
"corona" idea, which can be further tested in experiment. 

One direction is related with the predictions of the en- 
ergy /centrality dependence of the ridges and cones. Al- 
ready in the introduction we have shown in our sketch of 
the idea, depicted in Fig{TJ we indicated that since the 
geometry of the M-phase domain is quite different for dif- 
ferent collision energies. For RHIC the M-phase sits on 
the outside of the fireball (Fig. (a)), thus it get a maximal 
boost from the hydrodynamical expansion and may pro- 
duce ridges with rather narrow peaks in the azimuthal 
angle. For much lower energies - corresponding both to 
the SPS fixed target experiments and planned RHIC scan 
down) - one should find the M-phase only on the inner 
part of the fireball (Fig.(b)), which experience little hy- 
dro boost if at all. The logical prediction is then that 
no soft or hard ridges should be observed in this case. 
Standard modelling using the realistic geometry /density 
distributions should make those predictions quantitative. 

The "cone" is a different story: first of all, those do not 
rely on overall hydrodynamics and thus may well hap- 
pen at the very center of the firball and yet be observed. 
Second important distinction: "cones" are perturbations 
created by the "away-side" jet, and thus have completely 
different geometry/timimg. If the trigger jet is surface bi- 
ased, the away side jet has to fly through rather long path 
inside fireball, with its length varying roughly between its 
radius and diameter, 6-12 fm. As one can see from hydro 
solution we have shown in FigjT] , the M-phase starts at 
time zero at the edge and time about 5 fm/c at the cen- 
ter. Combining two observations together, one would see, 
that the away-side jet is travelling most of its long path 
in the M-phase. Furthermore, even at low collision en- 
ergies, when the M-phase occupies only the central part 
of the fireball, this conclusion is still fulfilled. The pre- 
diction then is that "cones" should not show very strong 
dependence on collision energy and centrality, in contrast 
to ridges. 

Let us now see how these ideas confront the available 
data. SPS experiments - CERES and NA49 - have not 
seen anything like ridges. The centrality and energy de- 
pendence of RHIC data we shave shown in FigjH] do in- 
deed suggest that this phenomenon is disappearing rather 
rapidly. On the other hand, both CERES and NA49 ob- 
serve away-side structures which are remarcably similar 
to RHIC data on the away-side, see review by Harald 
Appelshuser at the workshop [36] . We conclude that the 
picture in which ridges originate from spacial part of the 
"QGP corona", while cones are from its temporal part, 
at times 5-10 fm, is in qualitative agreement with the 
data. 

Finally, let me indicate one more direction of future 
studies: possible role of the unscreened electric fields in 
the M-phase in early-time hydro evolution. The pressure 
of the field, adding to (very low) pressure of the M-phase 
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may help to start hydro a bit earlier and help explain the 
HBT puzzle. Some studies of the kind (but with non- 
equilibrium fields rather than DMHD ones) have been 
made in Ref. [45] . 

IV. CONCLUSIONS 

In summary, we discussed two scenarios of the evolu- 
tion of extra energy /matter deposited by some fluctua- 
tions on top of the "Little bang" . In the scenario (A) we 
solved equations for propagation of sound with variable 
speed of sound in expanding matter, using the Hubble 
flow approximation. We have found that the rapid drop 
of the sound velocity in the "mixed" phase generates the 
secondary wave, which under certain conditions may be 
brighter and smaller in size, and thus is much better can- 
didate for the observed "cone" and "ridges". However 
this effect strongly depends on how sharply the speed of 
sound changes near T c , with the secondary wave being 
washed away if changes in the speed of sound gets are 
too smooth. 

In the scenario (B) we assumed that electric field re- 
mains unscreened for the duration of the M-phase (near 
T c ) and used dual Magnetohydrodynamics. We also 
found a potential for two expanding cylinders /cones. 
Furthermore, the speed of one mode can be zero which 
means that the flux tube of electric gauge field can be 
pressure- stabilized, by a (metastable) "coil" or current 
by magnetic charges in plasma. We have argued that 
this phenomenon is likely to stabilize the flux tubes in 
the near-Tc region, approximately at .STc < T < lATc. 
However this temperature interval does not account to all 
the time of the "Little Bang" at RHIC energies. Presum- 
ably the "coil" effect is still present at higher T, partially 
reducing the tube expansion. 

For the readers who may be surprised by similarity of 
the "double cones", let us remind that similar phenom- 
ena are well known in other fields of physics and may 
appear for multiple reason. In particular, in "dusty" 
strongly coupled electrodynamic plasmas double Mach 
cones have been experimentally observed (see e.g. [46 ): 
in this case these are sound and "shear" or hydroelastic 
modes which generate them. (Perhaps this option can 



lead to a "scenario C" not yet considered.) 

The very fact that we discuss two competing scenario 
should tell the reader that at the moment it is hard 
to tell whether they are robust enough to survive fur- 
ther scrutiny and explain pertinent observations. The 
"acoustical solution" is not quite robust, it does work 
only for rather sharp QCD transition which current lat- 
tice data do not support. Survival of electric field in 
the QGP corona is to be studied more, as well as pres- 
ence of metastable flux tubes. Clearly more theoreti- 
cal work is needed, including dedicated lattice studies of 
both monopoles and speed of sound. As far as experi- 
ment is concerned, it would be highly important to in- 
sure planning so that that jet correlations in question can 
be followed with sufficient accuracy during the planned 
RHIC scan toward the lower collision energies. Chang- 
ing the collision energy one changes the timing of the 
fireball eras, in a predictable manner, thus for the phe- 
nomena we proposed some quantitative predictions can 
be worked out and tested. 

Having said that, we emphasize that those theoretical 
and experimental studies seem to be very much justi- 
fied as both may potentially become important discover- 
ies. If the "acoustic scenario" (A) is the explanation, it 
would be a direct experimental signature of sharp quasi- 
lst-order QCD phase transition. If the DHMD scenario 
(B) would be confirmed, it would be quite significant find- 
ing, confirming reality of "QGP corona" in which physics 
is different from what it is inside, as much as it is on the 
Sun. It will be a big boost to "magnetic scenario" for the 
near-T c region J3TJ [32] . 
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